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EFFECT OF OXYGEN RECOMBINATION ON ONE-DIMENSIONAL 

AT HIGH MACH NUMBERS 



FLOW 


By Steve P. Heims 


SUMMARY 


A theoretical analysis of air flow in a channel in which oxygen 
dissociation and rec ombin ation occur has been made. The chang e] is viewed 
as a streamtube in the flow around a blunt body. The analysis is begun 
with the writing of the differential equation which gives the concentra- 
tion of atomic oxygen as a function of distance along the channel. The 
differential equation involves the reaction rate constant for the oxygen 
recombination reaction. This rate constant is evaluated theoretically 
from a formula due to Wigner, which yields a different result from simple 
collision theory. The equation for the atomic oxygen concentration thus 
obtained is solved together with the flow equations.- The equations may 
be solved by ordinary hand- computation procedures. An example is worked 
out to show the variation of the flow in a certain streamtube and its 
dependence on whether the oxygen reaction is “frozen/ 1 “in local equilib- 
rium," or proceeding at the finite rate indicated by the theory. The 
concept of a local relaxation length is employed. From inspection of 
the flow equations and the behavior of the cumulative lag of the dhe-ml csH 
reaction it is possible to judge without detailed numerical calculations 
whether changing one of the flow parameters brings the system closer to 
the “chemical equilibrium" or "frozen reaction" limit. 

An investigation iB made of the comparative relaxation times of 
the oxygen dissociation- recombination reaction in relation to molecular 
vibrations. A reason for interest in this is that it has -usually been 
assumed that vibrational relaxation occurs fast relative to chemical 
relaxation and therefore may be regarded as being in equilibrium. The 
present analysis indicates that it is not generally true that the vibra- 
tional relaxation times are short compared to the time characteristic of 
the chemical reaction. In this connection the generalization of the con- 
cept of chemical equilibrium constant is introduced for the case that the 
molecular vibrations are not in equilibrium. Some values of the relaxa- 
tion times are calculated and presented. A method is given to estimate 
the effect of vibrational lag when the vibrational relaxation times are 
relatively long. 
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* 

INTRODUCTION 


V- 

In the temperature range between 4-500° K and 7000° K and at densities 
lying somewhere between 0.1 normal atmospheric density and 10 times the 
normal density of the atmosphere, the oxygen of the air is largely in 

nascent form, while the nitrogen is mostly 
in diatomic form (see, e.g., ref. 1). This 
is the state of the air after it has come to 
chemical equilibrium behind a strong shock 
wave (Mach numbers 10 to 20, depending on 
the density of the cold air) . 

If subsequently this hot air flows around 
a blunt object, it will expand and cool and 
the oxygen atoms will recombine. The recombi- 
nation process is exothermic, giving off 5.12 
electron volts for every pair of oxygen atoms 
that recombine. Thus the reaction is a source 
of heat and will retard the coo lin g of the 
air. In this paper we study the effect of 
the oxygen recombination on the flow. Because 
of the difficulties involved in solving the 
flow around a blunt body even without chem- 
ical reaction, we focus attention on a single 
streamtube of the flow. Such a streamtube 
is approximated by the one-dimensional flow y 

in a channel of variable cross section. The 
cross section of the flow at first decreases, 
reaching a minim-urn which corresponds essen- 
tially to the sonic line, and then increases 
again (see sketch) . 

The content of the present study is indicated by the following 
outline: First a differential equation is derived, whose solution gives 

the fraction of the oxygen dissociated as a function of position in the 
channel. The formal solution is written down, but its numerical evalua- 
tion requires knowledge of the oxygen recombination rate constant and 
the value of the local pressures and temperatures. The rate constant is 
then discussed in the light of theory and experiment. Then the one- 
dimensional flow equations with chemical reaction and variable specific 
heat are given; a numerical method of solution of the flow equations is 
suggested. The three cases, zero reaction rate, infinite reaction rate, 
and the finite rate corresponding to the actual rate constant, are dis- 
tinguished a n d compared in a numerical example. Lastly the assumption 
of local vibrational equilibrium is exam in ed; some comparisons are made 
between vibrational relaxation times and chemical relaxation times and 
their practical and theoretical significance appraised. 
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SYMBOLS 


cross-sectional area 

cross-sectional area at throat of channel 
constant length in the Morse potential 

sum of the collision radii in collision between oxygen atom 
and third particle 

B,B x ,B(y) dimensionless parameter, defined by equation (12) 

C number of oxygen atoms, free or combined, per unit mass of 

fluid 

D dissociation energy of molecule 

E energy 

e base of natural logarithm 

F* statistical factor on which the recombination rate depends 

according to the Eyring theory 

f partition function 

g function related to enthalpy, defined by equation (28) 


A 

Amin 

a 


——s. fraction of states in oxygen-oxygen ground state to be counted 

jL®i as attractive 

h Planck* s constant 

H enthalpy per unit mass 

k Bolt7.ma.nn constant 

kj) dissociation rate constant 

kR recombination rate constant 

Ke equilibrium constant 

l length 

L length of channel 
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m 

mass of an atom; mass-flow rate through channel; also a 

# 


constant 


M 

mean molecular weight 


%t 

mean atomic weight 


k 2 

nitrogen 

r — 

NO 

nitric oxide 

- 

n Av 

Avogadro’s number 


*0 2 

number of particles of 0 2 per unit volume 

_ 

% 

number of particles of M per unit volume 


“0 

number of free oxygen atoms per unit volume 


0 

oxygen atom. 


0 2 

oxygen molecule 

.... * ' ' 

p 

pressure 

I 

p 

net number of 0 atoms liberated per sec per unit volume 


R o 

constant length in the formula for the Morse potential 

W 

r 

distance between two oxygen atoms 


S,S(y),Sx 

dimensionless measure, P — 

~ 

^0 \o(y) 


t 

time, sec 


T 

absolute temperature , °K 


u 

flow velocity 

~ . 

V 

relative velocity of two colliding oxygen atoms 


V(r) 

interaction potential between two oxygen atoms 


X 

fraction of oxygen atoms which are free, 110 

2no 2 +no 

■ 

*o 

initial value of x 

V 

xe 

value of x corresponding to local chemical equilibrium 

V 

y 

distance along the flow 
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% 

r(r) 

5 

q(T) 

A 


v 

P 

a 

T 


average vibrational energy per diatomic molecule 

efficiency of a triple collision in producing recombination 

ratio of change in vibrational energy to change in chemical 
energy, if vibrations are in equilibrium 

average efficiency of triple collisions in producing recom- 
bination 

relaxation length for chemical reaction defined by equa- 
tion (7) 

vibrational frequency of diatomic molecule 

density 

constant 

relaxation time of the chemical reaction 


Subscripts 


i 

at units of atmospheric pressure 

t 

o initial value 

si sea level 

v vibrational 

oo case of infinite rates of reaction 

(Thus Aoo is an approximate value of A where all the 
variables on the right-hand side of equation (8) are 
replaced by the infinite rate values.) 


Superscript 


o cold air 


EQUATION FOR THE FRACTION OF OXYGEN IN ATOMIC FORM 

As hot gas with partially dissociated oxygen enters a convergent- 
divergent channel, the cooling of the gas and the oxygen recombination 
reaction will proceed, causing a continuous variation of the fraction 
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of the oxygen which is dissociated. The purpose of this section is to 
derive and solve the equation expressing this variation. The flow treated 
is steady and one-dimensional. Catalytic reactions at surfaces are not 
considered. The solution obtained will depend on the value of the 
reaction rate constant, which will he treated separately in a later 
section. 


Derivation of Equation 



Since the flow is assumed steady, so 
that the condition at any one station is 
independent of time, we can write a time- 
independent conservation equation for 
free oxygen atoms. Let P equal the 
net number of oxygen atoms liberated per 
unit volume and unit time. Then 


PA£y 



dS 


( 1 ) 


simply states that the number of 0 atoms 
liberated per second in the volume element 
indicated in the sketch must equal the flux of free oxygen atoms across 
its surfaces. Here A is the cross-sectional area of the channel, 

dS an element of surface, u the local velocity of the fluid, and 
no the number of free oxygen atoms per unit volume. If the flow is 
one-dimensional, we have 



n-v u • dS = — (nouA)Zy 
v dy 


Defining x =no/(2no 2 +no) = mass ratio of the oxygen in free form to the 
oxygen in either the form 0 or 0 2 (mass fraction of oxygen in atomic 
foim) , we can write no =* Cpx, where C is a constant equal to the num- 
ber of oxygen atoms, either in the form 0 or 0 2 , per unit mass of fluid. 
Then equation (l) becomes: 


dx = _A_ ? _P_ 

dy mC puC 


( 2 ) 


where m = puA is the mass flow rate through the channel, a constant. 

Now P must be expressed in teims of the chemical rate constant 
and the concentration of the reactants. In order to dissociate, the 


t 


c 



NACA TN 4144 


7 


oxygen molecule must collide with another particle which is energetic 
enough to break up the oxygen molecule. The dissociation reaction can 
thus he written as: 


OsjfM 


»2 0+M 


where M is the "other" atom or molecule, which activates the oxygen. 
The rate of reaction will then he proportional to the number of colli - 
slons between 0 2 and M, so that one can write 1 

j = rate of dissociation = kjjfTjn^n^ (3) 

diss 

where k D is the dissociation rate constant, a function of temperature 
only. 


The inverse process, the recombination of two oxygen atoms, can only 
occur if some means is available to carry away the energy that the two 
separate atoms must lose to form a stable diatomic molecule. Two mech- 
anisms for carrying av?ay the energy are conceivable, a three-body colli- 
sion and radiation. The probability that this energy is given off by 
radiation is, however, negligibly small: Radiative (dipole) transition 

can occur only if one of the oxygen atoms is in the excited ^-D state, 
so the fraction of double collisions leading to combination is 

_ - 23.000 


10" 5 ^ e T Z. Here the 10 = factor is the usual ratio of duration 

<3 - 23,000 

R m 

of collision to the time required for radiation; the factor — e A 

, 9 

is the probability that one of the two atoms is in the D state; Z is 
a number much less than unity which gives the probability that the molec- 
ular state arising from this collision decays to the stable ground state 
of oxygen. The fraction of collisions where a third body is present is 
found from elementary kinetic theory to be of the order of 5x10 p a ^_; 
thus at atmospheric pressure and 5000° K a triple collision is 10 4 /Z 
more likely than a radiative collision. A general discussion of this 
point is given in reference 3, pages 400 to 402. 

In the case of a triple collision, the third particle takes up the 

kp 

excess energy and momentum. Thus the reaction is 0 + 0 + M — »0 2 + M, 
and 


Here the 10 


1 The implicit assumption is that the time between collisions (or more 
precisely, the correlation time) is long compared with the time it takes 
the 0 2 molecule to break up once it is activated. The latter time may 
be of the order of one period of vibration of the molecule or 2x10“ 14 

I rn -iq-IO 

sec, while the correlation time is larger than /— =— — sec. So we 

J 300 p o4 . 


sec. So we 


expect equation (3) to be valid for most cases of practical interest 
(see ref. 2) . 
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= rate of recombination ^ kp(T)nQ 2 nj4 (4) 

The form (4) of the rate law states the rate is proportional to the 
number of collisions between three bodies, two of which are oxygen atoms. 
As will be discussed in detail later (p. 17 ) , kp(T) has different values 
for different third bodies; thus when dealing with a gas mixture such as 
air, the right-hand side of equation (4) should be a sum of terms over 
different species (M x , M 2 , ...) in the gas. However at the present state 
of knowledge of the rate constant this appears to be an overrefinement, 
and the form (4) will be used. This form of the law has been verified 
experimentally for the recombination of like atoms in an inert gas atmos- 
phere at room temperature and near atmospheric pressure (refs. 4, 5, 
and 6) . 

According to all theories of reaction rates the rate constants k R 
and kp depend only on temperature, and do not depend on whether or not 
chemical equilibrium exists. They are related through the equilibrium 
constant Ke(T) = kp(T) /kp(T) . The equilibrium constant is a known 
function of temperature that can be calculated accurately from statis- 
tical mechanics; it permits us to eliminate either kp or kp from the 
equations. In the following analysis kp has been eliminated rather 
than kp . 

Since the net number of atoms liberated, P, is the balance between 
those freed by dissociation and the number of free atoms bound through 
recombination, one obtains with the help of equations (3) and (4) upon 
eliminating kp: 



P = k R n M (n 02 K e -n 0 2 ) 


(5) 


We can write equation (5) in terms of x, the fraction of oxygen atoms 
which are free, and Xe, the "local equilibrium value" of x. In other 
words, Xg is that value of x for which [no 2 K e (T) - uq 2 ] = 0. Then, 
noting that Kg = 2Cpx e 2 /(l-x e ) , one finds: 

P = -k R (T)C 2 p2n M + xj (x-xe) (6) 

Putting equation (6) into equation (2) gives 


dx _ -( x ~ x e) 
dy A 


with 
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1 + \ (, 

A u \ 1 -Xe / 

Equations ( 7 ) and ( 8 ) give the gradient of the quantity x at every 
point in the channel, and the integral of equation (7) gives the value 
of x at any point y. Equations (7) and ( 8 ) are valid no matter how 
far the deviation from chemical equilibrium and regardless of whether 
recombination or dissociation is the dominant process. 


Physical Interpretation and ’’Local Relaxation Length” 


The quantity A defined through equations ( 7 ) and ( 8 ) is called 
the "local relaxation length" of the reaction, and the time t = A/u 
is the corresponding relaxation time. 

In the special case that l/A and Xg are approximately constant 
throughout the flow region of interest, the integration of equation ( 7 ) 
yields an exponential decay to equilibrium: 


- Z 

x(y) - x e = [x( 0 )-xe]e * 


( 9 ) 


Here A is seen to be the distance required to reduce the deviation 
of x from its equilibrium value to one eth of what it is initially. 
Practically, A and Xe may be regarded as constant only when the density, 
temperature, and velocity vary slowly. In addition, A has some depend- 
ence on the deviation from local equilibrium, as can be seen by writing 
equation ( 8 ) in the form 


1 

A 


%(T) Cpn M 
u 



+ Xe) + (x-x e ) 


where it must be remembered that the number of third particles, ny, will 
usually have a linear dependence on (x-Xe) + xe. The dependence of l/A 
on (x-Xe) can always be safely neglected when (x-Xe) /x e is small compared 
to unity. 

Generally the physical situation is quite different from that 
described by equation ( 9 ) because the temperature and density may vary 
rapidly along the streamtube. Thus, for example, if the flow begins 
(at y = 0 ) in chemical equilibrium, but the parameters of the flow 
(p,T,u) vary continuously and rapidly along the streamtube, then the 
system, unable to keep up, will lag further and further behind as the 
flow proceeds. In this case, the deviation from equilibrium may increase 
rather than decrease with distance. Generally the tendency to lag behind 
equilibrium and the tendency to decay to local equilibrium exponentially 
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are operating simultaneously. Thus the qualitative "behavior of the 
system (whether it stays closer to the Infinite rate case or is nearer 
to the frozen case) is generally not determined "by the magnitude of the 
relaxation length alone, "but also by the magnitude of the gradient of x e 
and the gradient of A. These gradients can be expressed in terms of the 
gradients of temperature, density, and velocity. It is the relative mag- 
nitude of these gradients to the relaxation length that matters. "Whenever 
throughout the flow A« Le, where Le is the length within which Xe 
changes by a factor .of the order of one eth, the flow may be regarded as 
being in chemical equilibrium in the region of length L e . Deviations 
from equilibrium can be neglected in a channel of length L if every- 
where A« L, L e . On the other hand the flow is locally frozen when 
A/Le > > 1; it can be treated as frozen throughout a channel, if every- 
where A/L » 1. The more precise and complete criteria come from 
solving the differential equation for x. This will now be done. 


Solution of the Differential Equation 

Equation ( 7 ) together with equation .(8) expresses the conservation 
of oxygen atoms. In principle these equations should be solved simul- 
taneously with the other equations (such as mass and energy conservation) 
which the channel flow must satisfy. Although such a solution will be 
very tedious, it is always possible by numerical means "if the function 
k^(T) is known. 

However, an approximation to the solution is obtained by an itera- 
tion procedure which begins either with frozen flow or equilibrium flow. 
Here we shall for the sake of being definite deal only with the procedure 
beginning with chemical equilibrium. In this procedure one evaluates the 
factor l/A at the values p^, Too, Uoo and corresponding to the case 
of the reaction always being in equilibrium during the flow (infinite 
rate) . Similarly Xe can be replaced by x^" in equation (7) • The dis- 
tinction between Xg and x^, is a subtle one: x e is the equilibrium 

fraction of free oxygen atoms for the temperature distribution existing 
in the channel whatever the reaction rate; x^ is the special value 
of x e when the reaction rate is infinite. By this replacement it 
becomes unnecessary to solve equations (7) and (8) simultaneously with 
the other equations describing the flow, but one may first solve the flow 
for local chemical equilibrium, ignoring equation (7) , and then use the 
result to evaluate equation (8) and integrate equation (7). Let us call 
the approx ima te solution obtained thusly x( 1 ) . Then the error in x 
introduced by this approximation can be shown to be 

■ xti> ■ i (■* 1) s(y) 

n=o 


x 
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whe^e 


s(y) = (xg-O + X— [xU) -x^] 

\a 

is the error in first order. An upper limit to 5 may "be calculated ■by- 
evaluating xg and A at the temper at lore and density corresponding to 
chemically frozen flow. 

In the present approximation (l. e. , replacing x by xC 1 )) equation 
(7) is linear and the formal solution can be written immediately. This 
solution, which is the exact solution of the equation, 

dX^ 1 ) _ -[xO) -Xqq] 

d y 

is (see Appendix A) : 

x( x ) (y) =■ x^Cy) + e“ S ^) [xCoJ-x^fO) ] + e~ S ^ P e s *dx£, (10) 

^^(y) 


where the quantity S(y) is defined by S(y) = f — — , and S* 

^0 \o(y') 

the function of x^ which at every point is numerically equal to 
If one inverts the function x^(y) to obtain y(x co ) , then 
=S[y(x 00 )]. 

The last term in equation (10) vanishes when temperature and density 
(and therefore x^,) are constant along the flow. In this case it is seen 
that only the exponential decay term remains. On the other hand, if the 
fluid is in chemical equilibrium initially, the middle term vanishes 
leaving 

x( 1 )(y) = ^(y) + e“ S P e s *dxi (ll) 

^^(y) 

The solution (10) , which is constructed to be most accurate for near 
equilibrium flow, is seen to be exact also in the opposite limit of frozen 
flow (S = 0) . In practice 5 ‘and its derivatives are expected to be 
sufficiently small so that (10) will always be a fair approximation. 


is 

s(y). 
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Measure for Proximity to the M Frozen Flow” or 
"Equilibrium Flow" Limits 

The location of the solution relative to the two limiting cases 
("frozen" and "local equilibrium") is with the help of equation (10) 
given by: 

-S r r>3?oo(0) „ “I I 

-r-f j-r U(0)-x 0O (O)+ f e s djq,, (12) 

x(0)-x co (y) L J^y) -1 I 


Now B s 0 means the reaction occurs at an infinite rate, and B 1 
shows the reaction is frozen. It seems from equation (12) that one has 
to evaluate x(i)(y) completely in order to see whether B iB closer to 
zero or to unity. However, for some applications this is not necessary: 
If S(y) « 1, it follows from equation (12) that B(y) is near unity. 
This is a precise statement of the fact that if the relaxation length 
is long compared with the dimensions of interest, then the zero-rate 
approximation is valid. 

Another special application is to the flow leading from the region 
near the stagnation point of a blunt body around the edges to the sides 
of the body. In this case the reaction is expected to be in equilibrium 
at the initial point, taken in the stagnation regipn. Then it is clear 
from equation (12) that B has the property e“S = B = 1. Further B(y) 
can be given a simple geometric interpretation: It Is the ratio of the 

r x (°) a* 

two areas / e s dx and e b [x(0) -xn(y) ] . As Is to be expected from 

^Xco(y) 

the earlier discussion, the value of b depends on the shape of the 
function S^x^,) . 


B(y) - ^4^ 4 4 » 

x (o) -Xco(y) 
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In the Illustrative sketch above, the area of the rectangle enclosed 
by the heavy solid line represents eS*[x(0) , and the total area under 

one of the alternative dotted lines (which represent different functions 


S*) 


eS* 


r>Xc o(0) 
is / 

to B approaching e“ D and unity, respectively. For any convex shape 
such as 2, as well as for the straight line 3 } one has 1/2 < B < 1. In 
the last section use will be made of equation (12) to discuss the effect 
of body size and density on the solution. 


cLXoo. The extreme shapes 4 and 1 clearly correspond 

,-S 


The quantity B does not tell us whether both limiting cases are 
a good approximation or if they are both poor - it only tells of their 
relative merit. A more complete specification of the situation is 
obtained if the two parameters 


= x(y) -Xoo(y) ^ x(0) -x(y) 

Xco(y) * x(0) 


( 13 ) 


are employed. They represent the percent deviation from each of the 
limiting curves. They are related to B through: 


B 


\ + x (o) Ad ~ 

XooAoo J 


( 14 ) 


THE RECOMBINATION RATE CONSTANT 


This section is concerned with a detailed evaluation of knowledge 
concerning the recombination rate constant, a parameter upon which the 
flow depends. The reader who is not concerned with this aspect of the 
problem is invited to skip to page 21. 

At present there are no direct measurements available from which 
one can obtain the oxygen recombination rate k^(T) , for the temperatures 
of interest. However, considerable careful work (refs. 4, 5* and 6) has 
been done in measuring the rates of recombination of other diatomic gases 
(l 2 , Br 2 ) at room temperature in the presence of an inert gas. Recently 
a few results (refs. J, 8, and 9) have been published for I 2 and Br 2 
at temperatures between 1000° K and 2700° K in a shock tube. It was 
found that the recombination rate for these halogens is smaller at high 
temperatures than at low temperatures. 2 


2 After the calculations in the present paper were made, a shock- tube 
experiment using air was reported (ref. 10) . The author (S. Feldman) 
interprets the experiment to yield a lower limit for the recombination 
rates, which are about ten times larger than the rates given by equa- 
tion (21). 
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In order to interpolate or extrapolate from experimental results 
for one temperature to different temperatures or to a different diatomic 
gas or to different inert gases, a theory and some understanding of the 
mechanism of reaction is necessary. No altogether satisfactory theory 
for the oxygen recombination rate constant exists. There are, however, 
at least three methods available for estimating a rate constant from 
theoretical considerations. The three methods are (l) application of 
collision theory of chemical reactions, (2) application of the general 
theory of ’'absolute” reaction rates due to Eyring and co-workers (here- 
after simply referred to as the ’’Eyring theory”) , and (3) a method devel- 
oped by E. Wigner specifically for the combination reaction. Since each 
of the three theories can make a contribution to our understanding of the 
reaction, some comments are made about each approach; we finally use 
Wigner’ s theory to obtain a numerical value for the rate constant. Atten- 
tion will be focussed on the temperature dependence of the rate constant, 
because predictions based on the three theories give essentially different 
temperature dependencies if the usual approximations axe m ade in applying 
the theories. This is graphically Illustrated in figure 1. The graph 
will be further discussed at the end of this section. The disagreement 
between the theories may be due to the nature of the approximations made; 
it is possible that the rigorous calculation by Eyring theory, and further 
refinement of collision theory, would yield results in agreement with an 
exact calculation by Wigner’ s theory. 


Collision Theory 


The basic idea of collision theory (refs. 2, 6, and 11) is to count 
the number of recombinations by counting the number of triple collisions. 
Some arbitrariness comes in defining the effective radii of the colliding 
particles. This is resolved by approximating the atoms by hard spheres 
whose effective radii are taken to be the range of the interatomic forces. 
Experimental data such as viscosity measurements or spectroscopic meas- 
urements are often used to determine these radii (see ref. 12). Some- 
times experimental results are used to determine an efficiency factor, 
that is, the number, of reactions per triple collision. This adds another 
element of incompleteness to the theory. The advantage of collision 
theory is the simplicity and easy visualization of the concepts employed. 

If the duration of the collision is assumed to be inversely propor- 
tional to the relative velocity of the atoms before collision, then a 
rate constant of the form kft(T) = const «7 t is obtained. On the other 
hand, if this duration is assumed to be a function of only the inter- 
action of the two atoms, one obtains kj^(T) = const T. But no matter 
which way the duration is calculated, if the efficiency factor is regarded 
as a constant, the rate constant increases with temperature and disagrees 
with the iodine and bromine experiments. 

A serious shortcoming in the collision theory formulation lies in 
the neglect of the fact that the efficiency of a triple collision in 
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leading to chemical reaction is a function of temperature. Let tj(T) 
he this efficiency factor. Then k^(T) = const tj(T) ; T or 
kp(T) ss const t](T) *</T. It is of interest to get some idea as to the 
direction and the extent by which the introduction of ti(T) modifies the 
rate constant, so a very crude estimate of the efficiency of a triple 
col l ision is made in Appendix B. The efficiency factor is found to be 
the position average of f(r) where: 


r(r) 


V(r) 

kT 


A + * i |2gLf l 

\ kT 2 |_ kT J J 


( 15 ) 


Here V(r) is the value of the interaction potential between the two 
oxygen atoms at the point (r) where the triple collision is taking place. 
The dependence of r on l/T is very strong when V /kT << 1 (when r 
is close to zero), and becomes weaker as V/kT (and r) increases. The 
limiting foims are: 


r 


r = 


r a (r) 



i«l 

kT 


* 1 


»1 

kT 


(l6a) 

( 16 b) 


Formula (l6a) applies for large r where the potential between the two 
oxygen atoms is very weak, and (l6b) applies whenever the triple colli- 
sion occurs within the central part of the Morse curve. The value of 
t](T) is the appropriate average over the quantum mechanical states of 
the two oxygen atoms of: 


/ ^ r(r) - dr 

-TCP v(r) 

/ r 2 -$=£— 

» <0 v(r) 


(IT) 


In the averaging process, equation ( 17 ) , each volume element is 
regarded an equally likely site for a triple collision, except for a 
factor dr/v(r) which measures the time spent in the region dr during 
a collision. Note that for equation (l6a) the multiplicative factor 
r 2 /v is relatively large and for equation (l6b) it is relatively small. 
The different states to he summed over arise from different relative 
orientations of the angular momentum vector, which are possible between 
two atoms in their 3 P state. 
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V(r) 



Since t)(T) is an aver- 
age of r(r) , it will always 
decrease with, increasing 
temperature "because r(r) 
doeB. Just how rapidly it 
decreases will depend on 
the detailed shape of the 
potential. To illustrate 
these remarks consider the 
simple potential shown in 
the sketch; there Tx and r 2 
are independent of r. If 
V 2 /kT » 1, and Yx/kT « 1, 
then by the averaging pro- 
cedure of equation (17) 

ri(T) =( 1-a) rx+ar 2 rj ^ +a 



leading to 


k^(T) = const 



V T 


(18) 


Equation (l8) indicates that the efficiency factor could modify the 
simple collision theory sufficiently to change a rate constant with a 
positive square root temperature variation to one which decreases with 
increasing temperature, "but in any case the modification will be in that 
direction. 

A further small effect arises from the fact that at low temperatures 
the effective collision diameter on which the constant in equation (18) 
depends is not quite constant. It is shown in Appendix B that a factor 

proportional to T~^ is introduced into the recombination rate at low 
temperature on that account. Thus also this - factor tends to reduce the 
positive temperature coefficient which simple collision theory would pre- 
dict. Another approximation of collision theory (not further analyzed 
here) is to replace the duration of the collision by some average value, 
not taking into account the relative angular momentum of the two atoms. 
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A collision type theory which attempts to take these various effects 
into account properly would require an integration over two velocities, 
two impact parameters, and one spatial angle. The region of integration 
would again he given hy equation (Bl) . No such integration will he 
attempted here hut, as will he seen later on, Wigner’s theory does in 
effect just that, although using a language rather different from 
collision theory. 


Eyring Theory 


The approach used hy the Eyring theory (see refs. 13 , 14, and 15 ) 
is to diagram the total potential energy of the three atoms as a function 
of their relative positions. The diagram takes the form of an irregular 
surface similar to a system of mountains and valleys; the initial rela- 
tive position and velocity together with the interatomic forces cause 
the three atoms, represented hy a point on the diagram, to move in a cer- 
tain way. If the appropriate coordinate system is used, the motion of 
this point is identical with that of a mass point railing around under 
the action of gravity without friction in the corresponding mountain 
range. The combined state (0 2 ) is represented hy a certain valley, which 
is usually reached via a mountain pass. When the three atoms are in the 
particular relative position to each other, which is represented in the 
diagram hy a point in the neighborhood of this mountain pass, the group 
of atoms is referred to as an "activated complex"; before two oxygen atoms 
can combine, they must he part of an activated complex. 

According to the theory the rate constant is then seen as the pro- 
duct of a statistical factor F 4 " giving the concentration of the acti- 
vated complex, and the transmission coefficient k giving the probability 
that the activated complex will decay into the combined molecule plus 
third particle, k^(T) = kF 4 ". The activated complex is defined only for 
reactions requiring the addition of an activation energy; thus for two 
oxygen atoms which approach head-on, it cannot he defined because the 
potential energy curve has no maximum so that there is no activation 
energy required. However, generally two colliding oxygen atoms will have 
same relative angular momentum, and hence a centrifugal potential which 
yields a slight maximum in the potential energy curve. This fact permits 
the calculation of the recombination rate within the framework of the 
Eyring theory. We shall not attempt the complicated task of calculating 
the transmission coefficient for the oxygen recombination reaction. How- 
ever, some rather general observations can be made from Eyring’ s theory 
about the relative behavior of different third bodies: 

1. An inert third body is much less effective (k much smaller) 
in bringing about recombination than a third particle, such as an 
oxygen atom, which can form a chemical bond with one of the reacting 
atoms. 


2. The statistical factor F 4 " has a different temperature 
dependence for different third bodies. If the third body is a 
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hard sphere F* ~ */!?, the collision theory result; hut if there exists 
any - however slight - attractive force between the third atom, and 
the reacting atom/ then the temperature dependence lies between 
F*"~T and P f ~(T) 3/2 at the temperatures of aerodynamic interest. If 
the third body is a diatomic molecule, such as 0 2 or H 2 , then the 
temperature dependence lies between F + ~(T) 2 and F + ~ (T) 5 < 2 . The rea- 
son for these differences is that F* is proportional to the parti- 
tion function of the activated complex; in the hard sphere case the 
activated- complex is just a diatomic molecule, in the active atom 
case it is a triatomic molecule, and if the third body is diatomic, 
the complex consists of four atoms. The more atoms the complex 
involves, the larger the number of internal degrees of freedom and 
the sharper the dependence of the partition function on temperature. 

3. One cannot assert that the recombination rate constant will 
have a temperature dependence identical with that of F+ because k 
may vary with temperature. However, one does not expect the varia- 
tion of k with T to be so different for different third bodies as 
to wash out altogether the differences in temperature dependence due 
to differences in F*"(T), nor would it be expected that k decreases 
sufficiently rapidly as T increases to produce a rate constant with 
a negative temperature coefficient when the third body is a diatomic 
molecule. 


Wigner* s Theory 


Wigner (refs. 16 and 17) regards the three particles needed for a 
recombination as a member of a canonical ensemble (each member is three 
such point particles in a box, see ref. 18), and represents it as a 
point in (twelve-dimensional) phase space, where the 12 dimensions are 
the relative positions and momenta of the three particles. A certain 
region of this phase space represents the combined state, while the rest 
of the space represents the separated state. He then asks (and answers) 
the question: If initially the system is in ah uncombined state, what 

is the probability that after a time t it will be in the region repre- 
senting the combined state? He assumes that the particles follow the 
classical laws of mechanics and their interaction potential is a function 
of their relative position only , having no explicit time dependence 
(adiabatic hypothesis). These same assumptions also appear in the Eyring 
theory. For more detailed discussion of the assumptions underlying the 
theory and their validity see references 16, 17, .and 19 . For the case 
in which the third body is an inert atom, so that its interaction with 
the oxygen is that of a hard sphere, he obtains an especially simple 
upper limit for the rate constant: 


% - ^j ==z z [ dr(i«+2a 0M r 2 )V(r) 

pi s/flmkT J V< 


(19) 


0 
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where r is the distance between the oxygen atoms , ag^ is the sum of 
the collision radii for interaction between the oxygen and the third body, 
and V(r) is the interaction potential between the two oxygen atoms. The 
fraction g s /Sg.j_ is required to include the fact that the interaction 

between two oxygen atoms in their ground state depends on the relative 
orientation of their angular momenta. If several attractive states exist 
as is the case for oxygen (see, e. g. , ref. 3), the really correct proce- 
dure would be to add several terms of the form ( 19 ) , one for each attrac- 
tive state. But without considering the detailed shape of each of these 

potentials, one can estimate g B /Sg. on the basis of the formal Wigner- 

i ± 

Witmer rules (ref. 3). It Is noted that the values of g^ are constants 
and depend neither on temperature nor on the nature of the third body. 


The only dependence of equation (19) on the nature of the third body 
is through the ag^; the rate is larger, the bigger the inert atom. This 
is in quantitive agreement with experiments on I 2 and Br 2 in different 
rare gas atmospheres at room temperature. Further equation ( 19 ) predicts 

a temperature dependence for the upper limit k-^(T)— (T) 2 . That it is 
only an upper limit arises in the following way: In Wigner* s derivation 

of the rate constant, it is necessary to define the region of phase space 
corresponding to the combined state not only for the case that the third 
body is far away, but also while a] 3 three are in collision. Wigner 
shows that the correct definition of the latter region is the one which 
makes the rate constant a minimum, subject to the boundary condition that 
it reduces to the correct region when the third body is far away. Not 
having any simple procedure for obtaining the shape of this region, Wigner 
uses a simple and reasonable formula for the shape of the region; this 
will generally give too large a value for the rate constant. Moreover, 
the error introduced by using the wrong region increases with increasing 
temperature; consequently. It is expected that k^ decreases more rapidly 
as the temperature increases than equation ( 19 ) predicts (i.e., more 
rapidly than 1/ n/t) . 

For the present purpose, equation (19) will be used with sqm ® 2.60 A 
and a Morse potential function for V corresponding to the 3s g (ground) 
state of 0 2 : 


V(r) = D 


g(r-Rp) _ (r-Rp) ~~ 
e a _2e a 


( 20 ) 


with 


a = O.373 A , R 0 = 1.20 A , D = 5.12 e.v. = 59,000° K 

The fraction g s /^Si is estimated from the Wigner-Witmer rules to be 
l/lO. As was mentioned earlier, such an estimate involves some 
uncertainty. 
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The resulting rate constant is, if the third body is inert, 

10- (22i^)" a »ec- 


( 21 ) 


The corresponding relaxation times for various temperatures and densities, 
which will "be discussed in a later section, are shown in figure 8. 


Comparative Discussion of the Three Theories 


From collision theory a mechanical model or picture of the reaction 
is obtained, hut the theory is, in some respects, too incomplete, and in 
practice, requires the use of experimental data for application. 

From Eyring’ s theory some information is gained about relative rates 
with different third bodies. The difficulty in applying the Eyring theory 
to the present problem lies in the difficulty of calculating the trans- 
mission coefficient. The probability that the activated state leads to 
reaction is the probability that the particles end up not only in a cer- 
tain relative position (the two 0 atoms together), but also .with cer- 
tain relative momenta (the third particle carrying off most of the 
kinetic energy) . For this reason Eyring’ s use of a potential energy 
diagram in_ configuration space is not the most natural way to proceed, 
but rather the use of phase space is indicated where momenta and coordi- 
nates enter in symmetrical fashion. Also the use of the concept of 
activated complex seems not really Buited to the present problem, which 
involves practically no activation energy. 

Wigner, by working in phase space and avoiding the introduction of 
any activated complex, can express the rate constant in terms of an inte- 
gral, which he evaluates for the case that the third body is a hard 

— 1 

sphere. The resulting rate constant varies as T 1 2 ; the fact that the 
rate decreases as T increases is at least in qualitative agreement with 
the aforementioned I 2 and Br 2 experiments, although the da^a is best 

fitted by a T~s dependence. The modifications of this T 2 law which 
would result if the third body were a diatomic molecule rather than a 
hard sphere are not clear from Wigner' s integrals. The earlier discus- 
sion, in terms of Eyring' s activated complex, suggests that a less nega- 
tive or perhaps even a positive temperature coefficient would result. 

In the present analysis we shall make the simplest assumption, namely 
that all third bodies act as hard spheres, rather than assuming an 
activated complex. 

These arguments led the author to select for use in this analysis 
the theory of Wigner. The reasons for this choice are: 

1. Of the three methods his is the most naturally suited to 

the problem. 
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2 . Tie integrals obtained can actually be carried out, at least 
for a hard sphere third body. 

3. What experimental evidence there is, supports a negative 
temperature coefficient, and Wigner's rate constant gives a negative 
temperature coefficient, at least for a hard sphere third body. 

This concludes the critical evaluation of the three theories. It 
is of interest to look again at the graphical presentation (fig. 1 ) in 
the light of the above discussion because some more definite statements 
can be made: For oxygen in an inert gas presumably curve (I) is close 
to the correct form and (2) is a bad approximation. For oxygen in an 
atmosphere of N 2 the curve is expected to have a shape intermediate 
between (X) and ( 3 ) . The absolute value of curve (l) is obtained from 
equation (21) . No absolute values for the curves (2) and ( 3 ) are calcu- 
lated; they are arbitrarily taken to agree with (l) at 1000° K, permitting 
a comparison of the temperature dependence. These curves and their uncer- 
tainty serve to point up the need for further experimental measurements 
of the rate constants. 


FLOW WITH CHEMICAL REACTION 
The Equations 


Before one can integrate equation ( 7 ) the results (21) for the 
recombination rate constant must be inserted into equation ( 8 ) and A 
must be expressed in terms of x and A, rather than temperature (T) , 
density (p) , and velocity (u) . There is then obtained the following 
equation for x in terms of A, for any specified ch ann el shape A(y) : 

dx dy x-xg 
dA" ‘ dA A 

This equation is to be solved simultaneously with the r emaining flow 
equations to describe the flow completely. There are six variables: 

Xe, x, T, p, p, and u; hence, five additional equations are needed. These 
are given below. The equations allow for chemical reaction and variable 
specific heat but neglect transport phenomena: 


puA = m 

(mass conservation) 

(23) 

u du + ^ = 0 
P 

(momentum conservation) 

(24) 

u du + dH = 0 

(energy conservation) 

(25) 
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£ = (l+fx) 


~.x RT 

H 


(equation of state) 


K e (T) = 2Cp 


x e g 

l-x c 


(26) 

(27) 


The quantities M and f are defined by 


M = Sf iMi ; f = f q 2 + i f 0 


where 

^ number of moles of series i Qj 

total number of moles in cold air £nv° 

k K 

and _is the molecular weight of each species. For oxygen in an inert 

gas, M, f , and C are constant. However, for air both G and f depend 
on_the concentration of NO and thus are functions of temperature, but 
C/f is independent of the variable concentration of NO. The equilib- 
rium constant K e (T) and the enthalpy H are .computed from statistical 
mechanics; these functions have been tabulated by various authors (see 
refs. 1, 20, 21, 22, 23 , and 2k). 


Method of Solution 


Rather than attempting to solve all six equations ((22) to (27)) 
simultaneously, one can simplify the calculation by dealing with the 
frozen flow or equilibrium flow limit, depending on which is expected 
to be nearer to the case of interest, and then correcting for the finite 
rate case: 

For the frozen flow limit, one obtains the result by solving simul- 
taneously the four equations ((23) through (26)) using x = Xq = const 
in equation (26) . This differs from the well-known one -dimensional flow 
without reaction only in that care must be taken in evaluating the 
enthalpy H. The enthalpy at any point is obtained by adding to the 
initial enthalpy the change in enthalpy of each constituent due to the 
change of temperature, but not adding enthalpy change due to a heat of 
reaction, since no reaction is permitted. (See, e.g., the tables of 
refs. 21 and 22.) Then equations (27) and (22) can be solved in that 
order. Of the six variables, only x is sensitive to changes in the 
shape of the channel. If a higher approximation to the flow is desired, 
one then inserts the x(A) obtained from equation (22) into equation (26) 
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and, keeping p/p fixed, which is essentially equivalent to ignoring 
corrections to the enthalpy, obtains a new temperature at every point. 

In the other extreme case, equilibrium flow, one solves equations 
(23) to ( 27 ) simultaneously, replacing x in equation (2 6 ) by x^. The 
solution (ll) for (22) can then be employed . If a higher approximation 
to the flow is desired, one inserts the resulting x(a) - rather than x^ - 
into equation (2 6 ) and obtains a new T(A) . 

It is interesting to note that in both limiting cases, the properties 
of the flow (p, p, T, u) are a unique function of the cross-sectional 
area A, regardless of how slowly or rapidly the cross section varies 
along the channel. However, the corrections to the flow in the finite 
rate case are not unique functions of A, but through equation (22) 
depend on the channel shape, dA/dy; so for a finite reaction rate the 
flow is different for a channel which contracts (or expands) very sud- 
denly to a certain cross section than for one which contracts (expands) 
slowly to the same cross section. This dependence on the linear dis- 
tance arises because a characteristic length, the relaxation length, has 
entered into the problem. 

The simultaneous solution of equations ( 23 ) through (27) needed for 
the equilibrium flow is accomplished as follows: In analogy to the case 

of constant specific heat we introduce into equation (25) instead, of H, 
the variable g defined by: 3 

H « -S- E (28) 

g-1 P 

In general, g does not represent the ratio of specific heats. We then 
eliminate Xg from equations ( 26 ) and ( 27 ) and u from equations (24) 
and ( 25 ) . There are then obtained two equations involving only p, T, 
and p and one can find p(p) and T(p) from them numerically. The 
equations are 


Mp 

6RT 




0 = Ke(T) 


_P_ 

Po 


g- 1 'Po P 
(g“ 1 )o 


dp* 


(29) 


(30) 


3 This is the 7 * used by Gilmore, reference 1. He gives some numer- 
ical values for air. 
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Whenever g is constant along the path, equation (30) represents a 
polytrope : 



Usually g is a slowly varying function, and equation (30a) is a useful 
first approximation. The procedure used in the present calculation was 
to obtain the first p(p) from equation (30a), and to introduce this 
result into equation (29) to find T(p) . Then if it was found that 
g(p,T) did -indeed vary along the path, we fitted the slowly varying 
function g(p) "by a linear or quadratic function, thus 'avoiding numer- 
ical integration of equation (30) . This procedure converges very rapidly. 
Once p(p) and T(p) are known, one readily obtains without further iter- 
ation u(p) , p(A) , etc., from equations (23) and (24). In numerical 
evaluation the most inconvenient feature of this method is that in equa- 
tion (29) on the left-hand side one is subtracting two large numbers from 
each other to obtain a relatively small one and thereby losing some accu- 
racy. On the other hand a helpful feature is that the slopes of the 
right-hand side of equation (29), plotted as "a function of T, and that 
of the left-hand side are very different. 

As was noted in an earlier paragraph, it^ has up to now been assumed 
that the change in enthalpy dH in equation (25) is the same for flow 
with finite reaction as for the limiting case ( equilibrium or froz en ) 
on which the calculation is based. Now we show how the variation in dH 
may be evaluated. The enthalpy H (or the function g in eq. (30)) must 
be recalculated with the new x(A) and T(A) . This correction will have 
some (however small) effect on all the flow variables. In computing the 
enthalpy, it must be remembered that H depends on T, p, and x. This 
contrasts with the equilibrium situation, in which H depends only on T 
and p. To compute H for air with the help u3f available tables, it is 
useful to replace the actual history of the system by a different history, 
as is indicated in the accompanying sketch. This is permissible because H 


To find H or g at the 
point Ap, the history 
is taken to follow either 
the equilibrium curve or 
a (horizontal) line on 
which no reaction is 
taking place. 
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is a point function of T and x, and so it is .independent of the history 
of the gas expressed in terms of these three variables. With this modi- 
fication of the enthalpy determined, we repeat the solution and by this 
iteration obtain the correction to the flow variables. 

The enthalpy correction tends to be small not only because it is 

proportional to [x'^-x^,], but also because it is due to two partially 
cancelling effects: A change in T and a change in x in the opposite 

direction. These changes have opposite signs because p is kept fixed. 
The resulting corrections to the other flow variables for a given cor- 
rection to H are smaller the higher the temperature and flow velocity. 

The modification of the flow variables due to the correction of the 
enthalpy was found to be negligibly small in the example which follows. 


Example 


As a numerical example the following initial conditions were taken 
for the channel flow: 

— = 7*80 , 
p sl 


— = 0.336 , 
p sl 

where p g -| and p gl are the standard pressure and density at sea level 
and a, x is the rree-stream velocity of sound at an altitude of 25 km. 
These initial conditions correspond to the conditions of air having come 
to thermodynamic equilibrium behind a normal shock wave of Mach 14*7 at 
an altitude of 25 km* For the finite rate case, the length of the 
channel was taken to he 2*5 meters. 

Equations ( 26 ) and ( 27 ) were used with M = 28.8, f = O.OOO 765 , and 
C =s 0*22 n^ v /M^ . Thus the concentration of NO is neglected. For the 
rate constant equation (21) was used; the linearized version of the rate 
equation ( 7 ) was assumed so that the solution (10) could he employed* 

In figures 2 through 6 the results are plotted as a function of cross- 
sectional area. This area has to decrease sufficiently from its initial 
value so that the flow can reach sonic velocities; after that point the 
cross section increases again. 

It is seen from the graphs that for the frozen case (zero rate), 
the contraction of the streamtube at the throat is much less than for 
the infinite rate case. This shows that the oxygen reaction can change 
the shape of streamtubes in the flow. The effect reflects the relatively 


^ , 1.47 

a i 


To = 5300° K 



2 6 


KTACA TN 4144 


large equilibrium flow velocity at the throat compared to the slower 
frozen flow velocity. In the equilibrium case the enthalpy change between 
the initial point and the throat is comparatively large, being due to 
both thermal and chemical changes; the large enthalpy change requires 
a large increase in velocity, hence the greater contraction of the stream- 
tube in equilibrium flow (see sketch) . 



In addition, the temperature, 
density, pressure, and velocity are 
shown in'” the figures to be lower for 
frozen flow than for equilibrium 
flow when the flow is supersonic 
and when the area ratios are the 
same. Of course the variable x, 
the mass fraction of free oxygen, 
is always greater in the frozen 
case (xq) than in the equilibrium 
case (xoq) . 


The initial values of p, u, T, p, and x are taken to be the same 
for the zero rate and equilibrium case. From figures 3 and 4 it is seen 
that for both limiting cases the density falls^ below the free-stream 
value long before the pressure does. Although the curves for equilibrium 
and frozen flow appear to be very close to each other for both pressure 
and density, it must be remembered that their absolute values are small 
when A/Am -i -n is large and, consequently, the percent difference between 
the two cases may be considerable. In figures 2 through 6 all the curves 
have an infinite slope at the throat. This does not mean that the gra- 
dients are especially large in that neighborhood because dA/dy = 0 when 
A = Amin . and for any realistic streamtube dA/dy is very small in the 
neighborhood of A mi n. To compute x(A) (fig. 6) a particular shape had 
to be considered (see sketch) : 


A/Ajnin - 1-00 + 19-00 tanh 0.0586 (y/z ) 2 , Z = 75 cm (31) 



Cross section of channel, equation (31) 
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This shape has the property that for distances far from the throat 
(|y| >> Z) , the cross-sectional area approaches the constant value 
A = 20 A nrin » The value of B(y) turns out to he practically zero for 
small y and increases with increasing y. Its value for y = 230 cm 
is only 0.20. Thus in this example x(y) stays relatively close to 
x^(y), as is also seen in figure 6. The first-order correction to T 
is shown in figure 2. The enthalpy corrections (for T, p, p) discussed 
in the preceding section were found to he small, leading to practically 
no change from the infinite rate case. 

The comment made earlier about the lag of the reaction being cumula- 
tive can he examined from this example and, for this purpose, figure 7 
has been prepared. Even though the concentration of atomic oxygen lies 
relatively close to the infinite rate case, this figure shows that the 
percent deviation becomes progressively larger as the flow proceeds. 

Having considered a specific example, it is of interest to consider 
the effect of varying the parameters of the problem. One would intuitively 
expect that an increase of initial density would cause the reaction to 
stay still closer to chemical equilibrium (decrease the value of B) , 
because a higher density means more collisions and therefore faster reac- 
tions. Also one would expect that decreasing the length scale of the 
channel would cause the reaction to stay closer to the frozen limit 
(increasing the value B) . Both these expectations are borne out by the 
equations. (Qualitatively this can be seen from equations (8) and (12) 
alone; quantitatively density and distance must also be traced through 
the whole set of equations ( 23 ) to ( 27 ).) These facts show that one can 
simulate the chemical aspects of the flow at high altitude flight by 
using a small model at higher densities. 


WHEN CAN VIBRATIONAL RELAXATION BE IGNORED 


In the foregoing analysis the contribution of molecular vibrations 
to the internal energy of the gas has been included . However, it has 
been implicitly assumed that the equilibrium distribution of energy 
between vibrations and other degrees of freedom exists at all times. 

This assumption is investigated in the present section. Most of the 
results obtained apply to channel flow general l y; but some are applicable 
only to a channel flow where initially thermodynamic equilibrium exists; 
still others apply to the situation immediately behind a shock wave, 
before the gas has come to theimodynamic equilibrium. Whenever a result 
does not apply generally, the particular situation to which it does apply 
will be stated. 



28 


NACA TN 4144 


Two Alternative Conditions Permitting the Neglect of 
Vibrational Relaxation 


In order to find the circumstances under which vibrational equilib- 
rium can be assumed without introducing any significant errors , it is 
useful to define (see Bethe and Teller, ref. 21) a relaxation length for 
vibration, A v > which is the distance over which an initial deviation of 
the vibrational energy of the system from its equilibrium value will be 
reduced to one eth if the temperature of the gas does not change. This 
quantity is the exact analog to the chemical relaxation length A defined 
in equations (7) and (8) A gives the distance over which an initial 
deviation in x from its equilibrium value will be reduced to one eth, 
if the density and temperature of the gas remain constant. 

With this definition one can state in a quantitative way that vibra- 
tional equilibrium can always be assumed, if the vibrational relaxation 
happens sufficiently fast. Roughly, "sufficiently fast" means: 

^ A At 

— ^ <<1, <<1 everywhere (32) 

T A Ly 

Here Ly is the distance over which the equilibrium vibrational energy 
changes to one eth of its original value. If "in (32), much less than is 
replaced by much greater than, one has a sufficient condition permitting 
one to regard the vibrations as frozen. The more precise conditions, 
including cumulative effect, will not be discussed because for the pres- 
ent purpose only simple estimates based on (32) are required. These 
cumulative effects are analogous to those discussed in connection with 
chemical relaxation. The first of the conditions (32), t v /t << 1, states 
that the time scale of vibrational relaxation is so short compared to the 
chemical relaxation that, when initially both are out of equilibrium, the 
vibrations will already be adjusted before the chemical concentrations 
change significantly. The second statement says that vibrations will 
stay near equilibrium if the equilibrium position does not move away 
faster than the system* s rate of approach to equilibrium. 

Short relaxation times for vibrations provide a sufficient condition 
permitting the neglect of the effects of relaxation] but. it is not a 
necessary condition. This becomes clear if one considers low temperatures 
where vibrational relaxation times become very long] yet one can usually 
regard the flow as being in vibrational equilibrium simply because the 
fraction of the molecular energy tied up in vibrations is so minute, 
that whether vibrations relax rapidly or slowly will have no significant 
effect on the flow. Thus there exists a different, also sufficient but 
not necessary, condition allowing the neglect of vibrational relaxation 
which expresses the requirement that the amount of energy tied up in 
vibrations is sufficiently small: 
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AEv 

AH 


« 1 


(33) 


Here AEy is the change of vibrational energy for the case of vanishing 
vibrational relaxation time. The quantities AEV and AH are the changes 
of Ey- and H along a streamtube. When the condition (33) is satisfied, 
it is quite inconsequential -whether the vibrations are treated as frozen 
or as being in equilibrium. (The basis for the condition (33) is given 
in Appendix C.) 

In the remainder of this section these two conditions, (32) and (33 ) , 
are discussed in further detail; numerical estimates are made of the 
required quantities and applied to practical situations. 

The useful concept of ’’vibrational temperature" will be employed. 

It is defined in the following way: 

When there exists thermodynamic equilibrium, the vibrational energy 
per diatomic molecule is related to the temperature by 


€v( T ) 


hv 

e hv /'kT_]_ 


where h is Planck’s constant and v is the vibrational frequency of 
the molecule. In a nonequilibrium situation this relation between vibra- 
tional energy and temper at lire will not hold, but one can always find a 
number T v such that 


hv /kT v 
e -1 

It has been shown by Rubin, Shuler, and Montroll (refs. 25 and 26) that 
for all the cases of interest here, 4 the vibrational energy will be 
divided among the molecules according to a Boltzmann distribution, where 
T v plays the role of the temperature. Thus during the time when vibra- 
tional relaxation is occurring, the gas is correctly regarded as two 
systems (vibrations and the other degrees of freedom) in contact, each 
having a different temperature (T v and T, respectively) ; T v is called 
the vibrational temperature. 

4 More exactly, the situation considered by these authors is that of 
a system of harmonic oscillators which initially have a Boltzmann dis- 
tribution corresponding to a temperature T 0 , and which are contained in 
an excess of inert gas acting as a heat reservoir of temperature T. It 
is then proved that as time proceeds the system of oscillators continues 
to have a Boltzmann distribution, the temperature of which is the T v 
defined above. The coupling between the chemical reaction and vibrations 
may cause deviations from the Boltzmann distribution. 
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Numerical Values of Relaxation Times 


In the considerations concerning chemical relaxation behind shock 
waves, it was assumed by earlier authors (refs. 12 and 21) that the vibra- 
tional relaxation times are short compared to chemical relaxation times. 
Recent experimental evidence (refs. 10, 27, and 28) , however, seems to 
indicate that at high Mach numbers chemical equilibrium of oxygen is 
achieved before the vibrations have become active, so that condition (32) 
is not satisfied. 

Perhaps an intuitive reason for the earlier assumption is that an 
energy exchange between molecules sufficient to cause dissociation is 
much less likely than the small energy exchange required to produce vibra- 
tional excitation, and thus the time of relaxation must be longer. Here 
it must be remembered that the transition probability from 0 2 molecule 
plus M to two 0 atoms plus M is proportional to the density of final 
states, and this density is large for three free atoms, whose energies 
and momenta are not quantized. Furthermore, to achieve a vibrational 
state with energy nhv from the ground state, it requires n collisions, 
each of which must transfer a quantum hv of energy to the molecule in 
question; it requires more such collisions if the molecule loses some of 
its vibrational energy in between. However, "dissociation may occur in 
one step, although its probability is expected to be higher if the 
molecule is in a high vibrational state already (see refs. 27 and 29) • 

Numerical values for vibrational relaxation times were calculated 
by using the harmonic oscillator model (following Bethe and Teller) for 
vibrations, and the transition probabilities based on the rather satis- 
factory theory of Schwartz and Herzfeld (refs, 30 and 31) . Blackman’s 
experimental, data were used to estimate collision diameters (see ref. 

32). 


To make some numerical estimates of chemical relaxation times in 
air, equation (8) is used together with the recombination rate constant 
given by equation (21) . Above 6000° K it is found that' the chemical 
relaxation times are the same immediately behind the shock wave (x = 0, 
dissociation) and during the subsequent expansion process (x > x^, recom- 
bination) . At the lower temperatures there is some dependence on x. 

Thus for 5000° K and lower temperatures two values are given, the value 
for x = 0 and the one for x = Xq, the latter being put in parentheses. 
The chemical relaxation times are also expected to be different if the 
vibrations are frozen in instead of being in equilibrium. Both cases 
are tabulated in adjacent columns below and shown in figure 8. 
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T, 

°K 

p 

( t ^ 0 2 chemical reaction 

(tp) 0 2 vibrations 

^ t ^ ) ^N 2 vibrations 

Vibrations 
frozen at 
T v « 2300 ° K 

Vibrations 

active 

3,000 

10 

7-7 ( 3 . 7 ) 

8.1 ( 4 . 1 ) 

1.2 

40 

3,000 

.01 

330 ( 125 ) 

330 ( 125 ) 

1.2 

40 

4,000 

10 

.90 (. 45 ) 

1.38 (. 69 ) 

• 29 

16.5 

4,000 

.01 

12 (10) 

25 (20) 

• 29 

16.5 

5,000 

10 

.27 (. 18 ) 

.45 (- 30 ) 

.17 

5.6 

5,000 

.01 

.45 

1.25 

.17 

5.6 

7,000 

10 

.015 

.056 

.15 

2.2 

7,000 

.01 

.019 

.071 

.15 

2.2 

12,000 

! io 

.0012 

.0071 

.11 

• 57 

12,000 

.01 

.0012 

.0071 

.11 

• 57 


Note: The units for tp are microsecond atmospheres. The unit of p 

is normal atmospheric density. The quantity (tp) vibrations is i*ide- 
pendent of density or pressure , and (tp) chemical reaction s ^ 60 "becomes 
insensitive to changes in density and pressure at the higher tempera-* 
tures , 


The chemical relaxation time depends on the equilibrium constant 
through Xe and equation (27) - The frozen relaxation times were calcu- 
lated with the help of a generalization for vibrational nonequilibrium 
of the usual concept of an equilibrium constant of the 0 2 — > 2 0 reac- 
tion. By definition the equilibrium constant is a ratio of concentrations 
or the corresponding partition functions. When the distribution over 
vibrational states is that corresponding to a temperature T v rather 
than T, but that over other states is that corresponding to a tempera- 
ture T, the natural generalization fo r Ke(T) is: 

- Ke(T) ( 34 ) 

where f v (T) is the vibrational partition function corresponding to the 
temperature T : 

fv(T) =(l-e - hV/kT )' 1 

For vibrational equilibrium T v = T, while for frozen vibrations 
hv/kTy >>1 and f v (T v ) = For oxygen hv /k = 2300° K is the critical 
temperat lire, so that a vibrational temperature near ordinary room temper- 
ature leads to f v (T v ) = 1 . The evaluation of t for the frozen- 
vibrations case also requires that some assumption be made regarding the 
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dependence of tlie recombination rate constant on T v . The simplest, and 
not too unreasonable, assumption 5 is made that kp is independent of T v 
Wigner’s formula (21) is used throughout. 

In figure 8 the relaxation times are plotted as a function of tem- 
perature. The nearly horizontal regions into which vibrational relaxa- 
tion times fall are shaded in to point up how chemical relaxation times 
(based on Wigner’s theory) vary in relation to the vibrations. It is 
seen from this figure or the table that at temperatures in the neighbor- 
hood of 5000° K, the vibrational and chemical relaxation times are of 
the same order of magnitude, but that at higher temperatures the chemical 
relaxation times are relatively much shorter, and at lower temperatures 
they become longer than the vibrational relaxation times. In symbols; 

T / Tv ^ l , when T « 500Q° K 

t/tv «1 , when T > J000° K 

If the collision theory recombination rate used by Wood (ref. 12) 
is employed, one obtains still shorter chemical relaxation times. The 
region in which the two relaxation times are of the same order of magni- 
tude is of special interest: It suggests the possibility of studying 

the coupling of the two processes. This coupling between the two simul- 
taneous processes is sensitive to some details of the reaction mechanism, 
which otherwise would elude observation. For this reason further study 
of this coupling is considered to be of interest. 

The relative numerical values of r and r v have some implications 
for the flow immediately behind a strong shock, as well as for the flow 
of air around the comer of a rapidly moving blunt body. We consider 
the former case first. 

Behind shock waves the above result, that at higher temperatures 
chemical relaxation times are shorter than vibrational relaxation times, 
appears to be borne out by the unpublished results of Resler* s experi- 
ments with shock tubes (ref. 28) and also in line with corresponding 
experimental results for the dissociation of nitrogen obtained by the 
use of strong shock waves (refs. 27 and 29). Both the vibrational and 
chemical relaxation produce rapid cooling of the (translational temper- 
ature of the) gas and consequently a rapidly sinking equilibrium value 
for the degree of dissociation Xe and an increasing chemical relaxation 
time. As is suggested in the sketch, it is possible that x does not 
increase monatonically behind the shock wave but reaches a maximum and 
then decreases. This occurs as follows: At first the rate of approach 

to chemical equilibrium (eq. (7)) is rapid because it is proportional to 
(x-Xq) and because the temperature is high. As the reaction absorbs some 
of the thermal energy, the quantities x and Xe approach each other. 

5 The only such dependence on T v would arise if the third body in 
a triple collision were a molecule and its internal motion made it a 
better or worse energy absorber. 
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If the vibrational relaxation time 
is relatively long, the vibrations 
will absorb thermal (i. e. , random 
translational and rotational) 
energy, even after x and x e are 
very close to each other, and cause 
a continuing temperature decrease 
and decrease of xe* This permits 
the two curves x(y) and Xg(y) to 
intersect with a resulting maximum 
in x(y) . The experiments of 
Resler indicate such a maximum. 

Now in the expansion process, 
such as the channel flow considered here (see also ref. 8), the situation 
differs from the flow behind the shock in that vibrations are initially 
overexcited and the gas overdissociated, while the other extreme occurs 
immediately behind the shock. In the analysis of this expansion in the 
earlier part of this paper it was assumed that the vibrations are in 
equilibrium throughout the flow. Yet, since t v is not small compared 
to t, clearly the condition ( 32 ) is not satisfied. The justification 
for neglecting vibrational relaxation must then be based on the condi- 
tion (33) that the vibrational energy is sufficiently small. 


\ 



Application of the Sufficiently Small Vibrational Energy Condition 


In Appendix C it is shown that a practically equivalent way of 
stating the condition (33) , that the change of vibrational energy be 
small compared to the change in enthalpy, is: 


T) = 


| W + 


D 

5 Ax 


« 1 


(35) 


This condition states that the change of vibrational energy is smal 3 
compared to other changes in the enthalpy. 


For the special case of flow Immediately behind a shock wave , the 
vibrational energy is initially zero, so that Ae v = e v is the total 
energy per molecule received by the vibrational mode. From the analysis 
of Bethe and Teller it is known that the temperature behind a shock 
decreases, giving energy to the vibrations and the reaction. This 
exchange of energy would make q = 1, so that condition ( 35 ) is certainly 
not fulfilled. Actually q is usua ll y even larger than unity because 
the slowing down of the flow behind the shock front reduces the decrease 
in thermal energy (AT) which would otherwise occur. 
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It is quite different in the one-aimensional channel flow in -which 
one has vibrational equilibrium initially. In this case 
Ae v ~ kA T v = kAT . The equality pertains when there is no vibrational 
lag, so that AT V = AT. For oxygen D/5 k ~ 12,000° K; the condition (35) 
becomes 

i = 1 + 12,000 ^ » 1 

tj 2 ' AT 

In the particular example shown in figures 2 through 7, one sees from 
the graphs that Ax/AC ~ 1/5000° K; hence l/r\ 5 and the condition (35) 
is fairly well satisfied. 

A simple estimate of the enthalpy in this example indicates that 
the failure to consider finite vibrational relaxation times causes one 
to overestimate the changes in enthalpy by about 1 or 2 percent. From 
equation .(25) it then follows that the flow velocity changes a little 
more slowly than in the equilibrium case; the corresponding changes in 
the other flow variables can be easily obtained from equations (23) to 
(27) . The changes are so small that they could hardly show up on figures 
2 through 7» 


CONCLUDING REMARKS 


In conclusion we summarize the status of the problem of one- 
dimensional channel flow in the presence of recombination and dissociation. 

In principle, one only has to know the rate constant k R (T) , and 
one can, by numerical means, integrate forward the exact rate equation 
together with the flow equation, and all variables of interest will be 
given by the solution (see eqs. (22) to (27)). If another reaction is 
occurring simultaneously, another variable will be introduced, analogous 
to the mass fraction of oxygen in atomic form used for the dissociation- 
recombination reaction, and an additional equation must be considered. 

The simultaneous numerical evaluation of the flow and the rate of 
reaction is not necessary if the deviation from local chemical equilib- 
rium is not large. In this case one can compute the relaxation length 
(see eq. (8)) from the infinite rate approximation for the flow; then 
the solution for the mass fraction of oxygen dissociated takes on a 
relatively simple form, permitting numerical evaluation of the lag of 
atomic oxygen concentration behind equilibrium. This lag leads to 
readily computed modifications of the flow variables. One can express 
the amo unt of 1 n. g in terms of a parameter, whose numerical value identi- 
fies a flow which is essentially frozen or essentially in equilibrium. 

The chemical aspect of flow around a large body at low densities may be 
simulated by increasing density and decreasing size. 
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In particular it is noted that the shape of streamtubes, for given 
initial conditions, will he a function of the amount of lag behind chem- 
ical equilibrium. When the equations are applied to flow similar to that 
in a streamtube beginning near the stagnation region of a bluff body and 
leading around the body, it becomes apparent that the contraction of the 
streamtube in the sonic region is considerably greater when the reaction 
rate is very fast than when it is very slow. 

In experiments involving flow, the rate constant k-^(T) is not 
observed directly, but rather the relaxation length. This length depends 
not only on the temperature but also on local density, flow velocity, 
and extent of dissociation. Thus it requires a measurement of several 
quantities to determine the function kj^(T) at one temperature. Until 
further experimental data become available, however, Wigner* s theory of 
the recombination process is regarded as the best basis for estimating 
a rate constant. Collision theory fails to consider some important phys- 
ical features of the recombination process. The Eyring theory contributes 
information about the relative reaction rates for oxygen in an inert gas 
atmosphere compared to oxygen in an atmosphere of nitrogen molecules. 

It is not generally permissible to regard the vibrational relaxation 
as fast compared to chemical relaxation. In fact the present calcula- 
tions indicate that at high temperatures the rate of dissociation adjust- 
ment of oxygen is greater than the rate of vibrational adjustment. 

However, even when vibrational relaxation Is relatively slow, it may turn 
out to be entirely negligible because the energy involved in vibrations 
is small. The study of the coupling of the two relaxation times, when 
they are of the same order of magnitude, holds promise in yielding 
information about the reaction mechanism. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Sept. 10, 1957 
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APPENDIX A 

TEE SOLUTION OF EQUATION (7) 

To prove .that equation (10) in the text, 

x(y) = Xoo(y) + [x(0)-x 00 (y) ] + e“ S ( y ) P e s *dx£, 

^Xco(y) 

is the solution of the linearized form, of equation ( 7 ) : 

The standard form of the solution of the linear first-order equa- 
tion (7) is 


c(y) *5 P — e^( y * ^dy* + x(0) 

L Jq 


(Al) 


From the definition of S, one has 


r 5s eSCy'V . r y ^ SS. eSay. , f Uy V* S£ ** 
Jo \o Jo dy v Xoo(O) dx^. 


P T~ (^e 3 *)- 63 *] <kJo = e s x^ o -x OT (0)+ P 
Jx~(6\ J Jx 


2Qo(0) 


Uo) 


.s* 


^(y) 


Putting this result into equation (Al) gives equation (10) . 
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APPENDIX B 

COERECTIONS TO COLLISION THEORY 


The efficiency of a triple collision, r, is not, in general, a 
constant. A crude estimate of this quantity (ignoring the steric factor) 
might he made as follows: A triple collision is successful if the ini- 

tial velocities of the three particles are such as to lead to a final 
state in which the two oxygen atoms are hound: 

I (t) Native < V(r) (») 

The prime indicates final velocities. This condition on the velocities 
depends on r, the distance of the two oxygen atoms when the trip le 
collision is occurring. To simplify th i ngs we note yS - 

and replace (Bl) hy: 

imoVi ,2 <y(r) (B2) 

2 

The proposed model for the recombination process is then as follows: 
Oxygen atoms 1 and 3 are in collision (i.e., they are within a diameter 
of each other) . A foreign atom, 2, collides with 1 and reduces its 
velocity so as to satisfy (B2) . We ask what fraction of triple collisions 
satisfy (B2) . 

When it enters the potential field, particle 1 will increase its 
kinetic energy hy the amount V so the number which “while in collision 11 
have a velocity squared between v x 2 + 2V/m x and v x 2 + 2V /m x + dv x 2 is: 


f 


(f ^ + v ) 



m x v x 2 

2kT li dvn 2 
2 


At any one time the number in collision in a given velocity range is 
proportional to the duration of collision which varies as 

1/ vjv x 2 + 2V/m = l/v x . This approximation is evidently best when 
V/kT « 1. If the center of mass of 1 and 2 is assumed at rest then the 
relative frequency of triple collision of different velocities v^ is 

(mg/2jtkT) |v 3 {e-*^ Z 2 ^ s 0 finally we have for the average 

effectiveness of a triple collision: 
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r(r) S 


Vi ,2 < — 

mx 


mav 3 2 mxVi 2 " 

v 3 e 2 kT 2 kT dv^ dYx 2 


£ 


_ mev 3 2 _ mxYx 2 


v 3 e ^ 2kT dv 3 dvi 2 


all velocities 


For simplicity we take = m 2 » m, and then v 1 ,2 +v 3 ,2 = ^Vx 2 + + v 3 2 . 

To avoid the integration over angles we simply take 

. Then the integrations are elementary giving: 


V X * 2 = - 
2 


V IT + 


+ v 3 -= 


r(r) = 1 



(B3) 


When. even in the tail of the potential V/kT is greater than or of the 
order of unity, the effective size of the molecule decreases with rising 
temperature, because collisions taking place in the long Van der Waal 
tail of the potential curve will make an important contribution to the 
total number of successful reactions. To make some est imat e of this 
effect, we note from equation (B3) that r is a decreasing function 
as |v(r) j decreases. Now let us count only those triple collisions 
occurring within a radius such that on the average (over velocity) they 
have at least a probability r[V(r) /kT] = l/n of succeeding. For exam- 
ple, let us choose n such that in the region of interest, r can be 
approximated by some power law with an exponent s. Noting that the form 
of the potential for large r is v(r) = - c/r 6 , giving a maximum 


r 



then the diameter ~r mpv ~(e/kT) 1/6 . As the diameter enters into the 
expression for kp cubed, the effect of an increasing diameter is to 
introduce a factor s/c/kT. 
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APPENDIX C 

ANALYSIS OF CONDITION (33) FOR NEGLECTING VIBRATIONAL RELAXATION 


Generally one must look at all the flow equations, (22) to (27), 
to see the effect of vibrational disequilibrium. In a formal way these 
equations are unchanged, although an equation for the rate of change of 
vibrational energy (Ey) would have to be added. However the internal 
energy and thus the enthalpy H (eq. (25)) become a function of the 
amount of vibrational energy Ey and similarly the equilibrium con- 
stant Kg (eq. ( 27 )) will now depend both on T and Ey, and finally the 
chemical relaxation length A (eq. (22)) will depend on Ey. 

It is stated in the text that if the condition AEy/AH <<1 
(eq. (33)) Is satisfied, the effects of vibrational disequilibrium are 
unimportant. The quantity AEy is the change of vibrational energy 
assuming no vibrational lag, thus giving an upper limit to the change 
in Ey. We wish to show how this condition assures that neither H nor 
Ke nor A Ib importantly modified: 

That the modification of H is negligible is explicit in the con- 
dition (33) and needs no further discussion. The necessary modification 
of K e is expressed by equation (34) (i.e., Ke(T) is replaced by 
Kg(T, T v )). These two numbers are significantly different only when 
fy(T v ) and f v (T) are very different and, consequently, AEy is large 
and the inequality (33) cannot be fulfilled. (When one has a large 
temperature gradient, AH can become quite large, but A&y will increase 
proportionally. ) 

From equations ( 27 ) and (8) it is seen that the inverse chemical 
relaxation length is approximately linear in Ke and also depends on 
p, T, and u. The modifications required of p, T, and u to account 
for vibrational lag can be shown to be proportional to the modification 
in H. Thus if condition (33) is satisfied, neither Ke nor p, T, or u 
are much changed due to the finite vibrational relaxation time, and 
therefore l/X is not much changed either. 

Having pointed out that (33) is indeed a sufficient condition, we 
now show that condition (35) is practically equivalent to (33) : For a 

dissociating diatomic gas Ag of molecular weight of which a 

fraction X = n^ J /n^+2n^ e is in atomic form, the internal energy consists 
of random translational motion [Bp =3/2 RT/M^ 2 (l+ X) ] , rotational energy 
[E r = RT/Ma 2 (1- X) 3, vibrational energy [Ey = N/MA. 2 (l- *) e 3, energy due 
to dissociation [Eq = (ND/M^ e )X], and electronic energy E e q. One then 
obtains the following expression for the enthalpy per unit mass H by 
making use of the equation of state of the gas; 



4o 


WACA TN 4144 


H B E + S 
P 


Ep + Eg + Ey + Ejj + E e p_ + (l+X) 


RT 

Ma 2 


N 

%= 


+ + JEx + 

M Aa V 2 / 


E 


'el 


> (ci) 


The expression (Cl) does not assume either chemical or vibrational 
equilibrium. We make one simplifying assumption, namely, that changes 
in electronic excitation do not occur. This is certainly a good approx- 
imation at low temperatures, but the assumption must be dropped at very 
high temperatures. The characteristic temperatures (which appear in the 
Boltzmann factor) are 22,710° K and 27,610° K for atomic oxygen and 
nitrogen, respectively, while for molecular oxygen and nitrogen they are 
11,340° K and 71*580° K. For air, if the nitrogen is undissociated but 
a fraction x of the oxygen is dissociated, then 



Wow we have from equation (Cl) 


"ah 


1 

7 D 

■1 kAE + — Ax 

Ae v 


( 02 ) 


It is evident from equation (C2) that the condition AEy/AH << 1 is 


equivalent to 


Ae-y 


^ kAE + ^ Ak 
2 5 


« 1 . 


This is the condition (35) that was 


to be derived. 
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The constants for curves (5) and ® are 
chosen so that all three curves give 
the same value at 1,000° K . A 
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Figure 1.- Temperature dependence of the recombination rate by different 

theories. 
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Figure 2*- Temp e ratur e variation of air flowing through a converging-diverging channel with oxygen 

recombination. 
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Figure 3.- Density variation of air f loving through a converging-diverging cha n nel with oxygen 

recombination. 
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Figure 5.- Velocity variation of air flowing through a converging-diverging channel with oxygen 

re combine* ion * 
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Figure 6.- Degree of oxygen dissociation of air flowing through a converging-diverging channel with 

the sharpe given hy equation (31) . 
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Figure 7.- Lag of degree of oxygen dissociation behind its equilibrium value for air flowing through 
a converging-diverging nozzle with the shape given by equation (31) . 
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Figure 8 . - Relaxation times . 
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